One dimensional strongly interacting Luttinger liquid of lattice spinless fermions* 
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We consider the spinless fermion model with hard core 
repulsive potential between particles extended on a few lat- 
tice sites S. The Luttinger liquid behavior is studied for the 
different values of a hard core radius. We derive a critical ex- 
ponent Q of the one particle correlation function {c\cj) for an 
arbitrary electron density and coupling constant. Our results 
show that at the high density the behavior of fermions can 
be described as a strongly interacting Luttinger liquid with 
Q > 1. As a result, the residual Fermi surface disappears. 



The anisotropic one dimensional (ID) S=i Heisenberg 
model or the XXZ spin-chain is one of the simjjJe inte- 
grable models of the condensed matter physicall. The 
conception of Luttinger liquid introduced by HaldaneB is 
a very productive and useful for the description of one- di- 
mensional systemaj where power-law singularities in vari- 
ous correlators take place. The correlation lengths of spin 
correlation function and one-particle correlation function 
of spinless fermions at finite magnetic field and finite tem- 
perature have been calculated by means of a numerical 
study of the quantum transfer matrixQ. Authors have 
obtained a nonanalytic temperature dependence of corre- 
lation length of the spin correlator (uf (t|) and the Fermi 
momentum. They have found in this case a crossover 
temperature which depends on coupling constant and fill- 
ing. Another scenario of a behavior of the one-particle 
correlation function has been obtained in the framework 
of the integrable one dimensioual models with hard-core 
repulsive potentiaB . In Refs.El we discussed a strongly 
interacting Luttinger liquid state at high electron den- 
sity characterized by a large value of critical exponent 
for the momentum distribution function at T=0 ( > 1 
). In this case the residual Fermi surface disappears. A 
critical density ric (0(nc) = 1) separates a Luttinger liq- 
uid at a low electron density and an insulator phase that 
takes place at an extreme density nmax = jqr^ • Note, 
that a strongly interacting Luttinger liquid state and a 
high-Tc superconducting phase are realized at a small 
doping. 

In this Letter we consider the behavior of spinless 



fermions in the framework of the extended version of 
the XXZ spin-chain proposed by Alcaras and Bariev 
recentljti. As a concrete example we take the spinless 
fermion model with a hard core repulsive interaction. 
The model Hamiltonian contains a generalized projector 
that forbids configurations with two particles locating at 
distances less than or equal to 5. In this case the par- 
ticles interact only when they are at the closest possible 
positions (i,i+(5 +1, 5 is measured in units of the lat- 
tice spacing parameter). The case 5 — Q corresponds to 
the traditional XXZ spin-chain or the spinless fermion 
model. As we will see below, at high electron density 
the repulsive hard core potential leads to formation of 
a strongly interacting Luttinger liquid state without a 
residual Fermi surface. This example manifests a new 
many body state of ID fewnions, named as a strongly in- 
teracting Luttinger liquicu. As in Refs.cl, where this con- 
ception has been applied to different integrable models, 
we will discuss asymptotic behavior of the one-particle 
correlation function as a function of the electron density 
for several values of the coupling constant and the core 
radius. In the low temperature limit the critical expo- 
nents are evaluated of the scaling dimensions from finite 
size corrections to the energy spectra and can be calcu- 
lated using the thermodynamic Bethe equations. The 
behavior of XXZ-spin chain depends on the value of an 
anisotropic exchange interaction A: in the absence of an 
external magnetic field the point A = 1 separates the 
gap and gapless antiferromagnetic states. We shall con- 
sider both cases, namely A < 1 and A > 1, comparing 
a strongly interacting Luttinger liquid behavior in the 
insulator and metal states. 

The model can be considered as a generalization of the 
well-known spinless fermion model with the Hamiltonian 
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i=l 



^>5(4+iC» + c\ci+i -I- Anini+i+s)'P&, (1) 



* Dedicated to Renat Bariev's memory 



where c\ and Ci are spinless creation and annihilation op- 
erators of fermions at lattice site i, the hopping integral 
equal to unit and the coupling constant A is dimension- 
less, Vs is the projector forbidding two particles at dis- 
tances less than or equal to 5. By Ui ~ c[ci we denote 
the number operator for particles on site i. The system 
consists of N particles on the chain with L sites [L is as- 
sumed to be even). The Hamiltonian (1) is transformed 
into XXZ-chain using Jordan- Wigner representation . 
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Then, modftL (1) was exactly solved by the Bethe 
ansatz methoda. The two-particle scattering matrix is 
multiplied by an additional scattering phase shift due to 
the hard core potential 



S{ki, kj) = exp[-i(5(fci - kj)]SHiki, kj), 



(2) 



where Snikijkj) is the two-particle scattering matrix 
of the XXZ spin-chain 

g^g^^ fc ) = + '^^P^(fe» + fej) + Aexpifcj 
1 + expi{ki + kj) + Aexpiki 

At (5 = the two-particle scattering matrix (2) is re- 
duced to the one for the traditional XXZ chain. In the 
present Letter we consider two repulsive critical regimes 
0^ A < 1 and A > 1, introducing the parametriza- 
tion for the coupling constant A = cos 77 (0<77 ^ ^) and 
A = cosh rj , respectively. 

For 0^ A < 1 the charge rapidities Xj(j = 
1,...,N) related to the momenta of particles kj 

exp(ifcj) =~|~|~t[p~^) obtained by solving the 
Bethe ansatz equations 



( 



/ sinh i(Aj + irj) 
\sinh i(Aj — IT]) 
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sinh^(Aj- -Xi + 2ir]) 
.^^ sinhi(Aj - Ai -2i7?)' ^ ' 



where P — X^il^i the momentum. 

In terms of the rapidities A; the energy of the eigen- 
values is given by 



JV 



E = 2N cos 77-2 sin^ 77 



1 



cosh A, — cos 

j=i -I 



V 



(5) 



In the thermodynamic limit the ground state is de- 
scribed by the real rapidities {Xj}. Taking the ther- 
modynamic limit we obtain the integral equation of the 
Fredholm type for the distribution function p{X) for the 
variable A, 



p{X) + I dX'K<{X - A')p(A') = (1 - 5n)K<{X), (6) 

-A 



where the kernel K^{X) is, 



K<{X) 



1 



sin(i'77) 



2-K cosh A — cos(z^77) 



(7) 



The A-Fermi level denoted as A controls the band filling 
and the density of fermions is defined by 



dX p{X), 



(8) 



for ^ 71 ^ 77,0 and 



l-{l + 6)n= [ dXp{X), 



(9) 



for no < n ^ 7^,nax • The particle density Uq cor- 
responds to a 'half-filling' and rimax is an extreme 
density that corresponds to the full band. 

We remind the known results of the one-particle corre- 
lation function {c'{x)c{Q)) at T=0 obtained in the frame- 
work of conformal field theory. The correlation function 
shows an oscillatory behavior and power-like decay with 
a scaling dimension A' 



{c\x)c{0)) ~ cos(fcF2;)a; 



-2A' 



(10) 



here fc^? is the Fermi momentum. The momentum dis- 
tribution function close to /cf is determined by the expo- 
nent Q, 

{rik) — (nkp) — const\k — kp\^ sgn{k — kp), (11) 
where Q = 2A' — 1 and 



a V 2; 



(12) 



The long-distance power-law behavior of the spin cor- 
relator ((Tfcr|) is described by the critical exponent a; 
a=2C'^(A) is defined by the dressed charge C(A) at the A 
-Fermi level, C,{X) is the solution of the following integral 
equation 



C,{X)+ I dX'K<{X-X')C,{X')^{l-5n). (13) 

-A 



FIG. 1. Exponent a as a function of the electron density for 
rj = 0.1 (solid lines), 7r/4 (dashed); 7r/3 (dashed with points) 
and 5 = 0,1, 2, 3. Result for 5 = is plotted (dotted line) for 
comparision. 

FIG. 2. The exponent O , as in Fig. l.The line separates 
the strongly interacting Luttinger liquid state. 

The critical exponent a is analytically calculated for 

the densities n= (a = 2 ), tiq = 4^p^7iQ^ and 7T.niax 

{a = 2r7^ax) ■ This was done using Wiener-Hopf method. 
Solving numerically the integral equations (6), (13) and 
taking into account conditions (8), (9) we show on Figsl 
and 2 the critical exponents a and as a function of 
the density for some values of A (77 = 0, 1; 7r/4; 7r/3) and 
5 = 0,1,2,3 . The results of calculations obtained for 
the traditional XXZ chain (5 = 0) are plotted also for 
comparison in figures using dotted lines. In low density 
limit 71 ^ a converges to the value 2, which is the same 
as for non-interacting fermions. In the high density limit 
n — > rimax the value of a is independent on the coupling 
constant a = 2n'^^^ . In this case fermions are frozen and 
can be described by the second term in (1) . The model 
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Hamiltonian docs not depend on the coupling constant 
(so, we can choose it equals unit) and, as a result, the 
values of the critical exponents a and O are independent 
on A (see Figs 3 and 4 also). For 6 > a. density of holes 
changes in the interval from 6/{l+S) to 1 whereas the 
density of fermions varies from l/{l+d ) to 0. The hard 
core potential makes essential changes in the structure of 
the hole state. Note that this interaction also destroys 
the hole-particle symmetry. 

For 6 > the function is non-symmetric with respect 
to 'half-filling' point and the value of a is less than 1 
at n > 0.5. The correlation effects obtained due to the 
hard-core potential were most impressively displayed in 
the 8 behavior. According to (12) small values of the 
correlation exponent a leads to a large value of the 8 
correlation exponent (see Fig. 2). Comparing the behav- 
ior of 8 at (5 = and 6 ^ we observe an cixtrcinicly 
large value of 8 in the high density region which reaches 
^('^max ~ ^max)^ at maximal possible density rimax- In- 
stead of the trivial value of 8 = that takes place at 
(5 = we obtain i (|)^ at (5 = 1 , i (|)^ at (5 = 2 and 

^(^)^ at (5 = 3. The value of 8 increases with A and 
has a maximum value at a 'half-filling' for a large value of 

the coupling constant A. In Fig. 2 we separate by a dotted 
line the region of density that corresponds to a strongly 
interacting Luttinger liquid state with 8 > 1. The criti- 
cal exponent of the one-particle correlation function (10) 
also has a large value in this region ( 2 A' > 2). 

Using an analogous parametrization for the momenta 
exp(z/c,) — ^1" 2(^j+"?) obtained the Bethe equations 

for another region of the interaction coupling A = 
cosh > 1 in the form 



sin i {Xj + iff) 
sini(Aj- -iri) 
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= (-1)^ exp(-i^P) 

jj- sin^(Aj -A»-h2ir?) ^^^^ 



i=i 



sin i(Aj - Aj - 2ir]) 



For the distribution function p(A) and the dressed 
charge C(A) we have universal integral equations (6), (13) 
with a new kernel 



K>{X) = 



1 



sinh(j/77) 



2-77 cosh(i/ry) — cos(A) ' 



(15) 



FIG. 3. Exponent a as a function of the electron density 
for 77 = 0.1 (solid lines), 0.5 (dashed); 1 (dashed with points) 
and S = 0,1, 2, 3. Result for 5 = is plotted (dotted line) for 
comparision. 



FIG. 4. The exponent O , as in Fig. 3. The line separates 
the strongly interacting Luttinger liquid state. 



for different coupling constants A and S. The behavior 
of a is similar to the above case. Note that a minimal 
value of a is realized for an arbitrary coupling at a 'half- 
density'. A 'half filling' is a special point in one dimen- 
sional chains. As we note below, the gapless ground state 
for A < 1 transforms to the gap antiferromagnetic state 
for A > 1 . In Fig 4 the 6 - cusps characterize a gap 
phase for A > 1 , a residual peak takes place for A < 1 
at small 77 only (sec a solid line in Fig 2). This leads to 
the maximum of the value of 8 for arbitrary A and S (see 
the curves obtained for r] = 0.1; 0.5; 1 and ^ = 0, 1, 2, 3 in 
Fig. 4). 8 is an increasing function of the parameters of 
the interactions A and 6. According to numerical calcu- 
lation the critical density ric (8(nc) = 1 ) is less then no- 
The value of ric depends on the coupling constant A and 
the hard-core radius. 

A 'half-filling' point is a singular one in the behavior 
of the critical exponents. When A > 1 antiferromagnetic 
state with a gap takes place the 8 has a sharp maximum 
in this point. The width of the maximum is defined by a 
new scale, the gap Ag = 4sinhr7-^fc', where rj = ttK' /K, 
k and k' are modulus and complementary modulus of 
Jacobian elliptic functions and integrals, K = K{k) and 
K' = K{k') are complete and associated complete elliptic 
integrals of the first kind, respectively. The Ac trends to 
zero in the ry ^ limit. The elliptic integrals K and 
K' are connected with coupling constant i) and define a 
value of the gap. The gap in this state is a result of a 
large magnetic anisotropy. Thus, we can conclude that 
a strongly interacting Luttinger liquid is a state taking 
place near insulator or bad metal phases. It is possible 
also that a superconducting state has some features of a 
strongly interacting Luttinger liquid. 

In summary, using an exact solution of generalized 
XXZ-spin model we extended the calciilations of a 
strongly interacting Luttinger liquid state to the lattice 
spinless fermion system with a hard-core repulsive poten- 
tial. The system has two qualitatively different regions. 
In the low-density regime where the hard core potential 
is not essential we obtain a traditional Luttinger liquid 
. In opposite high density limit when the hard-core po- 
tential dominates the state is a strongly interacting Lut- 
tinger liquid with a large value of the critical exponent 
8 > 1. This behavior has been illustrated by the nu- 
merical calculations. It is shown, that the value of the 
critical density that separates these regions depends on 
both 5 and A. A main feature of this state is an absence 
of the residual Fermi surface. 

I.K. wishes to thank the support of the Visitor Pro- 
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In Figs 3,4 we show results of calculations of the cor- 
relation exponents as a function of the electron density 
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